A Multipole Based Treecode Using Spherical Harmonics for Potentials of the Form r-lambda

In this paper we describe an efficient algorithm for computing the potentials of the form r−λ where λ ≥ 1. This treecode algorithm uses spherical harmonics to compute multipole coefficients that are used to evaluate these potentials. The key idea in this algorithm is the use of Gegenbauer polynomials to represent r−λ in a manner analogous to the use of Legendre polynomials for the expansion of the Coulomb potential r−1. We exploit the relationship between Gegenbauer and Legendre polynomials to come up with a natural generalization of the multipole expansion theorem used in the classical fast multipole algorithm [2]. This theorem is used with a hierarchical scheme to compute the potentials. The resulting algorithm has known error bounds and can be easily implemented with modification to the existing fast multipole algorithm. The complexity of the algorithm is O(pN logN) and has several advantages over the existing Cartesian coordinates based expansion schemes.